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Ever  increasing  acceptance  of  electric  vehicles  relies  on  better  operation  and  control  of  large  battery 
packs.  The  individual  modules  in  large  battery  packs  cannot  have  identical  characteristics  and  may 
degrade  differently  due  to  manufacturing  variability  and  other  factors.  Degraded  battery  modules  waste 
more  power,  affecting  the  performance  and  economy  for  the  whole  battery  pack.  Also,  such  impact  varies 
with  different  trip  patterns.  It  is  beneficial  if  the  performance  gain  can  be  evaluated  online  prior  to 
replacing  certain  battery  module  regarding  to  actual  driving.  In  this  paper,  the  research  objectives  are 
two-fold.  First,  we  propose  an  on-line  battery  module  degradation  diagnostic  scheme  using  the  intrinsic 
signals  of  battery  module  equalization.  Without  need  for  additional  sensors  or  offline  tests,  this  scheme  is 
cost-effective  for  constructing  and  updating  the  battery  pack  “health  map”  in  real  time  during  the  vehicle 
operation.  Secondly,  based  on  the  derived  battery  health  map,  the  Worthiness  of  Replacement  (WOR)  for 
certain  modules/cells  is  proposed  to  evaluate  the  performance  of  the  battery  pack  for  customer  specified 
trip.  Such  evaluation  index  provides  a  quantitative  measure  for  module  replacement  and  battery  pack 
swapping  Simulation  study  is  performed  to  validate  the  proposed  ideas  using  an  actual  driving  cycle 
recorded  for  a  commuting  trip. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electric  vehicles  (EV)  require  large  battery  packs  with  high 
energy  and  power  densities  to  become  a  competitive  choice  of 
transport.  There  has  been  extensive  research  on  battery  equaliza¬ 
tion  and  management  systems  for  EV  [1],  Various  battery 
management  systems  have  been  developed  to  achieve  cost- 
effective  and  fast  voltage  equalization,  simultaneously  with 
battery  state-of-charge  (SOC)  and  health  monitoring  [2,3],  The 
inter-module  voltage  and  charge  equalization  is  critical  for  healthy 
and  efficient  operation  of  large  battery  pack  [4],  Battery  equaliza¬ 
tion  can  relieve  the  negative  impact  of  degraded  modules,  although 
not  reversing  the  degradation  [5],  However,  with  degraded 
module(s)  present,  there  is  significant  power  loss  during  the 
equalization  process  due  to  the  increased  internal  resistance,  and 
such  loss  is  also  affected  by  the  driving  cycle  (or  the  trip  pattern).  It 
is  of  EV  customers’  interest  to  evaluate  the  (negative)  consequences 
of  degraded  module(s)  and  decide  if  replacement  is  necessary.  The 
decision  of  replacement  thus  relies  on  a  combination  of  both 
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diagnosis  of  physical  conditions  and  economical  payback.  More 
specifically,  one  desires  to  know:  1 )  how  much  energy  loss  would 
be  induced  for  a  specific  driving  cycle  (trip),  based  on  the  knowl¬ 
edge  of  module  characteristics  within  the  battery  pack,  and  2)  for 
the  given  battery  pack,  what  energy  efficiency  it  would  be  by 
replacing  certain  degraded  module  or  cell  with  a  healthy  one. 
Recently,  we  have  investigated  these  two  issues  considering  the 
internal  resistance  as  the  primary  indicator  for  battery  health  [6]. 

The  solution  to  the  first  problem  can  be  facilitated  via  online 
acquisition  of  battery  health  map,  and  more  specifically,  the 
module-wise  internal  resistance  map  [7],  In  this  study,  we  propose 
to  identify  the  internal  resistances  of  individual  modules  via  the 
signals  during  module  equalization  process.  From  the  system 
identification  standpoint,  this  approach  takes  the  driving  cycle 
power  demand  as  input  excitation  to  identify  the  parameters  (i.e. 
internal  resistances)  needed  for  battery  health  mapping.  Without 
need  for  any  extra  sensor,  this  method  has  the  cost-effective 
advantage. 

For  the  second  issue,  this  study  proposes  a  quantitative  measure 
for  the  evaluation  of  module  based  replacement,  named  as  the 
“Worthiness  of  Replacement”  (WOR).  For  a  given  driving  cycle,  the 
WOR  of  a  specific  battery  pack  is  defined  as  the  ratio  of  the  SOC 
change  of  the  current  battery  pack  to  that  after  replacing  certain 
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module  with  a  healthy  module  of  nominal  specifications.  Any 
performance  evaluation  for  vehicle  powertrain  relies  on  the  use  of 
driving  cycle.  Such  evaluation  can  be  performed  for  standard  (e.g. 
DOE  or  EPA)  driving  cycles  or  any  user  specified  trip  (e.g.  from 
onboard  GPS)  because  the  driving  cycle  pattern  can  strongly  affect 
the  power  demand  profile,  and  in  turn  the  current  flow  through  the 
whole  battery  pack  and  the  equalization  currents  among  the  indi¬ 
vidual  modules  [8],  Recently,  there  have  been  significant  efforts  in 
using  trip  information  and  modeling  to  optimize  the  energy 
management  for  HEV  (hybrid  electric  vehicles),  PHEV  (plug-in 
hybrid  electric  vehicles)  and  EV  [9],  In  this  study,  we  will  evaluate 
the  WOR  concept  via  actual  driving  cycle  data. 

The  above  two  aspects  of  work  are  validated  with  an  EV 
model,  parameterized  like  a  light-duty  delivery  vehicle.  An 
example  service  trip  during  the  morning  rush  hours  is  used  for 
illustration,  which  consists  of  both  local  road  and  freeway  road 
segments.  An  existing  battery  equalization  strategy  is  imple¬ 
mented  [10],  and  a  subspace  identification  method  is  applied  to 
identify  the  battery  module  internal  resistances  via  the  charging 
and  discharging  signals.  The  battery  internal  resistance  map  thus 
obtained  is  also  SOC  dependent,  which  makes  the  WOR  evalua¬ 
tion  more  accurate. 

The  remainder  of  this  paper  is  organized  as  follows.  The  EV  and 
battery  models  are  presented  in  next  section.  Section  3  describes 
the  battery  module  equalization  scheme  adopted  in  this  study. 
Then,  the  subspace  identification  based  internal— resistance  esti¬ 
mation  is  presented  in  Section  4.  Simulation  results  are  presented 
in  Section  5,  and  the  work  is  concluded  in  Section  6. 

2.  Electric  vehicle  and  battery  model 

2.1.  Electric  vehicle  configuration  &  modeling 

The  EV  model  used  for  this  study  is  assumed  to  be  driven  by  a  4- 
pole,  560  V,  permanent  magnet  synchronous  motor  (PMSM)  with 
a  gearbox  transmission.  The  rated  speed  is  5000  rpm.  There  are  two 
different  battery  pack  configuration  which  will  be  described  in 
Section  5.  The  torque  limits  were  set  to  140  N-m  during  driving  and 
160  N-m  during  braking,  respectively.  In  order  to  derive  the  motor 
driving  torque  and  motor  current,  the  vehicle  propulsion  dynamics 
is  needed.  Based  on  the  propulsional  dynamics,  the  wheel  torque  is 
[11] 


I'wh  =  (Fad+Frr  +  FacC)r  (1) 

where,  Fa d  =  l/2pACav2  is  the  aerodynamic  resistance  force, 
face  =  Ma  is  the  acceleration  force,  and  Frr  =  p.„Mg  is  the  rolling 
resistance  force.  M  is  the  vehicle  mass  and  a  is  the  acceleration. 
is  the  rolling  resistance  coefficient,  p  is  the  air  density,  A  is  the 
effective  frontal  area,  Cd  is  the  drag  coefficient,  and  v  is  the  vehicle 
speed.  The  driving  torque  at  the  motor  shaft  is 

=  <2> 

where,  Tm  is  the  motor  shaft  torque,  TWh  is  the  torque  at  wheels,  r  is 
the  wheel  radius,  G  is  the  gear  ratio,  and  r)s  is  the  gear  efficiency.  In 
this  study,  the  vehicle  parameters  are  set  as  follows:  M  =  1560  kg, 
G  =  11,  Cd  =  0.19,  p.rr  =  0.0048, 77g  =  0.95  and  r  =  0.3  m. 

2.2.  Battery  models 

It  is  assumed  that  battery  is  the  only  power  source  for  the 
simulated  EV  system,  and  the  lithium-ion  battery  is  assumed.  There 
are  two  battery  models  used  in  this  study:  one  is  the  simulation 


model  based  on  detailed  battery  physics,  and  the  other  is  a  simpli¬ 
fied  model  for  subspace  system  identification  based  internal 
resistance  estimation. 

To  develop  the  virtual  plant  of  electric  vehicle  propulsion  for 
simulation  study,  a  detailed  lithium-ion  battery  model  is  adopted 
from  MATLAB  Simulink®  SimPowerSystem™,  which  consists  of 
nonlinear  time-domain  models  for  charging  and  discharging 
processes,  respectively  [12,13]. 

For  discharging  process  (i*  >  0), 


E  =  Eo-K-Q^i-KQ^it+Aexpi-Bit)  (3a) 

For  charging  process  (i*  <  0), 


E (3b) 

where  £  is  no-load  voltage,  Eo  is  constant  voltage,  K  is  polarization 
constant,  i *  is  the  low-frequency  current,  i  is  the  battery  current,  it  is 
the  extracted  capacity,  Q  is  maximum  battery  capacity,  A  is  expo¬ 
nential  voltage  and  B  is  exponential  capacity. 

For  the  online  estimation  of  battery  module  internal  resistance, 
a  simplified  linearized  R—C  battery  model  is  used.  Among  the 
various  battery  models  we  adopted  a  battery  model  from  [14]  for 
internal  resistance  estimation,  as  shown  in  Fig.  1.  In  the  R—C  model, 
Cs  and  Cl  are  surface  capacitance  and  bulk  capacitance,  respectively. 
Rs,  Rl  and  Rq  are  surface  resistance,  bulk  resistance  and  series 
resistance,  respectively.  In  general,  Rq  and  Rs  are  much  smaller  than 
bulk  resistance  RL.  Similarly,  surface  capacitance  Cs  is  much  smaller 
than  bulk  capacitor  Cl. 

In  order  to  compute  the  SOC  trajectory,  we  have  adopted  the 
following  dynamics  [15] 


SOC(k+l)  =  SOC(k) 


Voc  -  y/v, o2c  -  4(Rint)Tm  •  Wm  •  Vmssn^ 
2(Rint) -Qb 


(4) 


where,  Voc  is  the  open  circuit  voltage  of  the  battery,  Tm  is  the 
electric  motor  torque,  wm  is  the  electric  motor  speed,  ijm  is  electric 
motor  efficiency,  Rjnt  is  the  internal  resistance  of  the  battery  which 
is  equivalent  to  Rl  in  the  battery  model  in  Fig.  1,  Qb  is  the  battery 
capacity  and  SOC  (k)  is  the  state  of  charge  at  the  k-th  time  instant. 
The  SOC  calculated  with  coulomb  counting  is  the  charging  and 
discharging  of  the  cells  during  a  specified  amount  of  time.  Eq.  (3) 
reveals  that  the  SOC  change  depends  on  the  operation  mode  and 
load  demand.  As  the  motor  speed  and  motor  torque  increases  the 
battery  SOC  drain  increases.  Typically,  the  healthy  upper  limits 
for  battery  SOC  is  80-90%,  while  the  healthy  lower  limits  is  about 
20—30%.  In  this  study,  we  have  made  a  relatively  conservative 
choice  on  the  battery  SOC  range  from  80%  to  30%  [16,17], 


Fig.  1.  Battery  equivalent  circuit. 
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3.  Battery  equalization 

3.1.  Battery  equalization  scheme 

As  the  voltage  in  a  single  cell  is  quite  low,  for  high  voltage 
applications  like  EV,  battery  modules  and  packs  are  made  via  serial 
and  parallel  combinations  of  the  cells  [18].  Variations  are  inevitable 
in  the  internal  resistances  of  the  battery  cells  and  modules.  The 
variation  in  the  internal  resistance  is  caused  by  extrinsic  or  intrinsic 
cell  properties,  contact  resistance  among  cells,  or  temperature 
gradient  because  of  improper  thermal  management  accounts  for 
battery  imbalance.  The  variation  in  the  internal  resistance  causes 
variations  in  the  SOC  and  degrades  the  cell  [19].  The  battery  life  can 
be  significantly  affected  by  the  imbalance  among  modules  or 
cells  due  to  the  associated  over-charging  or  over-discharging. 
Hence,  equalizing  the  modules  or  cells  is  important  for  acquiring 
the  maximum  achievable  power  and  extending  the  battery  life 
[20,21], 

In  this  study,  the  voltage  equalization  scheme  by  Lee  and  Cheng 
is  adopted  [10],  with  a  two-module  situation  shown  in  Fig.  2.The 
voltage  of  each  module  determines  the  direction  of  energy  transfer 
between  the  two  modules  with  proper  operation  of  the  MOSFET 
switches  Qi  and  Q2.  Li  and  I2  are  two  uncoupled  inductors,  while  C\ 
is  an  energy  transferring  capacitor. Vbi  and  Vfe  are  battery  voltages 
for  modules  (or  cells)  1  and  2,  respectively.  For  normal  condition, 
Vci  =  Vbi  +  Vb2.  The  ideal  condition  is  Vbi  =  Vb2,  although  this  can 
seldom  be  the  case  in  realistic  operation. 

With  assistance  of  Fig.  3,  two  operational  scenarios  are 
described  in  the  following.  For  the  PWM  operation,  T$  is  the  period 
and  D  denotes  the  duty  cycle. 

1)  VBi  >  VB2:  For  duration  DTS,  Qi  is  turned  on,  and  capacitor 
Ci  transfers  energy  to  Vb2  during  the  same  period  inductor  Li 
stores  energy.  For  duration  (1  -  DTS ),  Qi  is  turned  off,  D2  turns 
on,  and  the  capacitor  energy  is  supplied  to  Module  #1  and  L2 
charged  by  Vb2.  In  the  switching  duty  cycle  (Ts),  as  shown  in 
Fig. 3,  T,  -  To- DTS,  while  T2  -  Ti  =  (1  -  DTS). The  operation  can 
be  described  by  the  following  equations. 

a)  For  t  e  [To,  Ti],  Qi  is  turned  on,  which  implies 
VB1(t)  =  tid!^  with  iu(T)0  =  Jo  (5) 


VB2(t)  =  ~L2^f  +  ^  J  !L2«dT  With  iL2(To)  =  70  (6) 

To 

b)  For  t  £  [7i,  T2],  Qi  is  turned  off  and  Q2  is  turned  on,  which 
implies 


cSi  cL1T:t! 

2=r^=r 
Tn  r 

rv - 7- 

j  ■ 

Q2  L  i=- 

"i  E2 

Fig.  3.  Typical  switching  waveform  of  the  battery  equalizer  when  VBi  >  VB2. 


).i,«2+2./illWd, 


with  iu(Tt)  =  IP 


(7) 


VB2(t)  =  -L2d!^  with  il2(T2)  =  fP  (8) 

Eqs.  (4)  through  (7)  describe  voltages  VBi  and  VB2  dynamics 
during  one  duty  cycle  of  PWM  operation.  Ip  is  the  peak  current 
during  the  operation  cycle  DTS.  The  currents  in  the  inductor  7li  and 
7l2  are  given  by 

fu  =  [5  (j^D2  +  '4lj~'/B1(l  -  D)2)]  Ts  (9) 


=  [5  (Vc1L2Vb2<1  -  d)2+^0  -  £>)2)]  Ts  (10) 

By  varying  the  switching  frequency,  the  equalization  scheme 
can  be  implemented  in  continuous  current  or  discontinuous 
inductor  current  mode.  During  the  equalization  process,  the  weak 
cell  is  charged  by  a  strong  cell  to  balance  the  energy  level  of  both 
the  cells.  In  our  study  battery  1  is  at  higher  energy  level  than 
battery  2. 

2)  Vbi  <  Vb2:  The  operation  is  controlled  by  the  switching  Q2  and 
D 1  in  similar  fashion  as  Case  1,  and  similar  equations  would 
apply. 


Fig.  2.  Battery  equalization  cir 
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3.2.  Worthiness  of  replacement 

The  concept  of  worthiness  of  replacement  (WOR)  is  applied  to 
the  battery  pack  after  the  battery  is  equalized  and  the  trip  is 
complete.  As  mentioned  earlier  though  equalization  is  unavoidable 
but  WOR  projects  the  exact  loss  of  energy  in  the  trip.  This  analysis 
would  be  the  decision  making  factor  for  projection  of  completion  of 
the  trip  as  the  onboard  energy  is  calculated.  The  WOR  is  defined  as 


is  linear  in  B  and  D.  The  predictor  is  formed  from  past  inputs,  and 
Eq.  (14)  is  an  output  error  model.  Therefore,  if  the  system  operates 
in  open  loop,  one  can  consistently  estimate  B  and  D. 

y(t|B,  D,  xQ)  m  c(ql-  a)  x08(t)  +  c(ql-  a)  Bu(t)  +  Du(t) 

(16) 


WOR 


SOC  Change  of  Current  Battery  Pack 
[Battery  Pack  SOC  Change  with  Certain  Module(s)  Replaced] 


(11) 


Eq.  (10)  can  be  enhanced  by  some  averaging  operation,  for 
example  for  actual  or  predicted  trips  of  certain  number  of  days. 
Considering  the  stochastic  factors  in  traffic  conditions  and  driving 
behaviors,  such  averaging  can  help  to  obtain  more  reliable  result  for 
evolution  purpose. 


Eq.  (15)  is  linear  in  Xo  =  x  (0)  and  <5(t)  is  the  unit  pulse  at  time 
t=  0. 

Matrices  A  and  C  are  estimated  from  extended  observability 
matrix  for  r  observable  rows,  i.e. 


4.  Subspace  estimation  of  internal  resistance  MAP  of  EV 
battery  pack 

In  order  to  obtain  the  module-wise  internal  resistance  map  of 
battery  pack  for  WOR  evaluation  and  online  monitoring,  we 
propose  to  estimate  the  internal  resistance  with  the  battery 
equalization  signals,  i.e.  the  associated  voltage  and  current 
measurements,  during  the  actual  trip  (i.e.  with  the  driving  cycle 
input). 

In  this  study,  the  subspace  identification  [22]  method  is  applied 
for  the  relevant  parameter  estimation  based  on  the  current  and 
voltage  measurements  in  the  equalization  circuit.  Based  on  a  state- 
space  battery  model,  subspace  parameter  estimation  permits  esti¬ 
mation  of  static  battery  parameters  based  on  input— output  data 
(fin  and  V0).  The  process  can  be  carried  out  online  during  the  vehicle 
operation  so  that  actual  change  of  battery  characteristics  can  be 
reflected  [23], 

The  subspace  identification  theory  is  briefly  reviewed  in  the 
following.  For  the  m-input,  p-output,  n-th  order  system  with 
discrete-time  model 

x(t+l)  =  Ax{t)  +  Bu{t)  +  w(t)  (12a) 

y(t)  =  Cx(t)+Du(t)  +  v(t)  (12b) 

where  x(t )  are  states,  u(t)  and  y(t)  are  input  and  output, 
respectively,  w(t)  is  the  process  noise  or  disturbance,  and  v(t)  is  the 
measurement  noise.  System  (11)  can  be  diagonalized  via  similarity 
transformation  as 

x(t  +  1)  -  T^ATx(t)  +  T~*Bu(t)  +  w(t)  (13a) 

y(t)  =  CTx(t)  +  Du(t)  +  v(t)  (13b) 


where  T  is  invertible  and  defined  by 

m  =  TMt)  (14) 

The  output  can  be  estimated  by  assuming  that  A  and  C  are  fixed,  i.e. 

y(t\B,D)  =  c(qI-A)  Bu(t)  +  Du(t)  (15) 


G  has  dimensions  (pr)  x  n.  C  can  be  estimated  by 

C  =  G(l:p,l:n)  (18) 

Similarly,  A  is  estimated  by 

C(p  +  1  :  pr,  1  :  n)  =  G(1  :  p(r  -  1),  1  :  n)A  (19) 

After  estimating  the  initial  state,  the  prediction  error  minimi¬ 
zation  method  can  be  used  to  determine  the  unknown  parameters 
in  the  state-space  model.  In  this  study,  the  target  parameter  is  the 
internal  resistances  for  the  respective  batteries.  The  estimation 
quality  can  be  evaluated  with  the  prediction  error 
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Motor  speed  during  driving. 


Motor  torque  profile  during  driving. 

Fig.  5.  Driving  cycle  and  torque  profile  for  the  simulation  example. 


z(t,d*)  =  y(t)  —y(t\6*)  (20) 

fori  <  t  <  N  the  error  sequence  is  filtered  by 

sF(t,d)  =  L(q)s(t,d)  (21) 

where,  L(q )  =  1  in  this  case.  For  parameter  estimation,  the  cost 
function  is  defined  as 
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where,  /(■)  is  a  quadratic  function,  ZN  =  [y(l),u(l),y(2),u(2),  -, 
y(N),u(N)]  is  the  vector  of  the  input— output  data  pairs.  The  least- 
square  method  is  applied  for  parameter  estimation: 


(24a) 


(24b) 


dN  =  6, y(zN)  =  argmin  VN(d,ZN^  (23) 

where,  JeDMc  RN,  and  DM  is  the  set  of  parameter  vectors.  This 
procedure  of  estimating  the  required  parameter  8  is  prediction- 
error  identification  method  [22], 

In  order  to  apply  the  foregoing  prediction-error  identification 
method  to  battery  equalization  process,  the  associated  state-space 
model  is  obtained  as  follows.  Based  on  the  two  operational  modes 
described  in  Section  3.1,  the  averaged  state-space  model  was 
derived,  as  presented  in  the  Appendix.  When  battery  1  is  at  higher 
energy  level,  battery  1  charges  battery  2.  Referred  to  Fig.  2,  the 
averaged  state-space  model  for  battery  1  is 


where,  Xi  is  the  voltage  across  capacitor  Csi,  x2  is  the  voltage  across 
capacitor  Cli,  X3  is  the  current  through  inductor  Li,  and  X4  is  the 
voltage  across  capacitor  Co-  Similarly,  the  state-space  model  for 
battery  2  is: 


Table  1 

Initial  conditions  of  two-module  simulation  study. 

Case  no.  Battery  no.  Voltage  (V)  Capacity  (Ah)  SOC  (%)  Internal 

_ resistance  (0) 

1  B-l  100.8  20  80  0.075 

B-2  100.8  20  80  0.075 

2  B-l  100.8  20  80  0.075 

B-2  100.8  20  65  0.30 


298 


R.  Ugle  et  al.  /  Journal  of  Power  Sources  223  (2013)  293-305 


Time  (sec) 

Fig.  6.  Battery  module  SOC  trajectories  for  with  Case  #1  with  the  example  driving  cycle. 


Fig.  7.  Two-module  battery  SOC  for  example  driving  cycle  for  Case  2. 
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Road  in  Thiensville,  while  the  destination  is  3200  North  Cramer 
Street  in  Milwaukee. 

The  proposed  idea  is  evaluated  with  simulation  study  based  on 
the  EV  model  described  in  Section  2.  The  Matlab,  Simulink  and 
SimPowerSystem  are  used  for  the  simulation.  We  used  a  lithium- 
ion  battery  model  as  a  power  source  for  the  electric  vehicle 
[12,13],  For  the  PMSM  motor,  the  driving  torque  can  be  obtained  by 
applying  the  driving  cycle  (speed  profile)  to  the  propulsion 
dynamics.  The  torque  limits  were  set  to  140  and  -160  N  m.  The 
current  limits  for  the  battery  side  were  set  from  250  A  to  -250  A  on 
the  battery  side.  The  driving  cycle  had  both  long  and  short 
segments  of  speeds  [25],  The  torque  and  speed  profiles  for  the 
actual  driving  cycle  are  given  in  Fig.  5,  which  will  later  be  used  to 
simulate  two  and  three  module  batteries  with  different  cases.  The 
concept  of  WOR  in  two  different  scenarios:  a  two-module  battery 
pack  first  and  then  a  three-module  battery  pack. 


where,  the  states  are  same  as  those  in  Eq.  (23). 

After  obtaining  the  values  for  the  parameters  from  the  subspace 
state-space  estimation  we  calculate  the  values  for  each  of  the 
resistances.  Rest.  Rcli  and  Rq\  are  identified  internal  resistances  for 
battery  1  and  Rcs2.  Ren,  and  R02  are  internal  resistances  for  battery 
2.  Thus  the  internal  resistances  of  individual  modules  can  be  esti¬ 
mated  throughout  a  specific  trip. 

5.  Simulation  results 

Simulation  study  has  been  conducted  to  evaluate  the  proposed 
method.  The  driving  cycle  adopted  is  the  one  used  in  an  earlier 
study  [24],  as  shown  in  Fig.  4.  The  trip  origin  is  124  West  Freistadt 


Table  2 

Initial  conditions  of  individual  battery  modules  for  simulation  study. 


Case  no.  Battery  no. 


B-l 

B-2 

B3 

B-l 

B-2 

B-3 


B-2 

B-3 


Voltage  (V)  Capacity  (Ah)  SOC  (%)  Internal 

resistance  (O) 


64.8  20 

64.8  20 

64.8  20 

64.8  20 

64.8  20 

64.8  20 

64.8  20 

64.8  20 

64.8  20 


80  0.054 

80  0.054 

80  0.054 

80  0.054 

72.5  0.081 

65  0.108 

80  0.054 

80  0.081 

80  0.108 
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Fig.  8.  SOC  trajectory  for  Case  1  of  three-module  battery  pack. 
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Fig.  9.  Battery  SOC  trajectories  for  three  modules  for  Case  2  of  Table  2. 


5.1.  Simulation  study  for  two-module  battery  pack 

Two  battery  modules  are  assumed  serially  connected.  For  each 
module,  the  capacity  rating  is  20  Ah,  respectively.  The  upper  and 


lower  limits  of  battery  SOC  are  set  to  be  80%  and  30%,  respectively. 
Each  module  consists  of  4  sub-modules  which  have  7  cells  in  series, 
i.e.  28  cells  in  total,  with  nominal  voltage  3.6  V  [26].  The  voltage 
across  the  battery  pack  is  201.6  V  [27],  i.e.  100.8  V  for  each  module. 


Fig.  10.  Battery  SOC  trajectories  for  three  modules  for  Case  2  for  intermediate  equalization. 
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Two  cases  are  simulated  with  the  settings  described  below  (also 
summarized  in  Table  1,  with  B-l  and  B-2  representing  battery 
modules  1  and  2,  respectively): 

Case  1  Both  modules  are  in  the  nominal  condition,  i.e.  having 
internal  resistance  of  0.075  Q.  Starting  from  80%  SOC,  both 
modules  follow  the  same  discharging  trajectory  as  shown 
in  Fig.  6.The  final  SOC  is  53.52%. 

Case  2  Battery  1  is  in  nominal  condition,  while  battery  2  is 
degraded  with  internal  resistance  of  0.30  O  [28],  As  battery 
2  is  degraded,  we  assume  the  starting  SOC  level  is  65%. 
Such  difference  in  SOC  level  results  in  the  equalization 
process  described  earlier.  In  Fig.  7,  we  observe  that  the  two 
modules  equalize  during  the  first  645  s  of  the  driving 
cycles.  The  energy  from  B-l  charges  B-2  to  increase  its  SOC 
level.  The  two  modules  are  equalized  to  61.8%  SOC,  and 
then  both  modules  discharged  in  identical  fashion  down  to 
33.45%  SOC  by  the  end  of  the  driving  cycle.  Though  the 
driving  cycle  is  complete  in  Case  2  there  is  a  considerable 
amount  of  energy  loss  from  B-l  supplied  to  B-2  so  as  to 
equalize  which  could  have  been  used  for  driving.  For  Case  2 
in  particular,  it  can  be  observed  that,  due  to  the  relatively 
large  internal  resistance  of  B-2,  there  is  significant  amount 
of  energy  loss  during  the  equalization  from  B-l  to  B-2.  The 
battery  pack  SOC  changes  by  26.48%  in  Case  1  (both 
modules  healthy),  while  by  39.05%  in  Case  2.  Such  differ¬ 
ence  leads  to  the  WOR  of  1.47  for  B-2. 


5.2.  Simulation  study  for  three-module  battery  pack 

Simulation  is  then  performed  for  a  battery  pack  with  three 
modules  in  series.  The  capacity  for  each  module  is  20  Ah.  The 
battery  pack  is  composed  of  three  modules.  Each  module  consists 
of  two  serially  connected  sub-modules  of  nine  cells  each,  i.e.  18  cells 
in  series.  This  configuration  results  is  64.8  V  across  each  module. 
The  overall  voltage  across  the  battery  pack  is  194.4  V.  The  internal 
resistance  for  ideal  module  is  assumed  to  be  0.054  Q. 

Three  cases  are  simulated,  with  settings  listed  in  Table  2,  with  B-l, 
B-2  and  B-3  standing  for  battery  modules  1, 2  and  3,  respectively.  For 
Case  1,  all  the  modules  are  in  nominal  condition,  i.e.  with  internal 
resistance  of 0.075  Q  and  80%  initial  SOC.  For  Case  2,  B-l  is  in  nominal 
condition,  B-2  has  internal  resistance  of  0.15  Q  with  72.5%  initial  SOC, 
and  B-3  has  internal  resistance  of  0.3  Q  with  65%  initial  SOC.  This 
case  is  intended  to  make  B-2  and  B-3  more  degraded  modules.  Case  3 
is  similar  to  case  2  except  that  the  starting  SOC  is  identically  at  80%. 

For  Case  1,  all  the  modules  demonstrate  identical  discharging 
trajectories  as  shown  in  Fig.  8,  with  the  final  SOC  of 46.26%.  In  Case  2, 
as  shown  in  Table  2,  B-l,  B-2,  and  B-3  have  different  initial  SOC  levels, 
which  results  in  equalization  among  the  three  modules  from  the 
start  till  1005  s,  as  shown  in  Fig.  9.  In  particular,  B-l  has  higher  SOC 
than  B-2,  and  B-2  has  higher  SOC  than  B-3.  The  equalization  of  these 
three  battery  modules  follows  the  same  logic  as  described  earlier,  i.e. 
the  equalization  decision  is  made  for  every  pair  of  neighboring 
modules  based  on  the  comparison  of  the  relevant  voltages,  which  is 
done  by  triggering  the  associated  switches.  In  the  scenario  assumed 
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for  Case  2,  current  is  discharged  from  B-l  to  B-2,  and  current  is 
discharged  from  B-2  into  B-3.  All  the  modules  equalize  their  SOC  to 
45.57%.  Afterwards,  the  three  modules  jointly  discharge  to  35.34% 
SOC  by  the  end  of  the  driving  cycle.  Although,  in  Case  2,  the  driving 
cycle  can  be  completed  with  the  onboard  battery  power,  there  is 
a  considerable  amount  of  SOC  from  B-l  supplied  to  B-2,  and  that 
from  B-2  supplied  to  B-3.  Further,  more,  due  to  the  relatively  higher 
internal  resistances  of  B-2  and  B-3,  there  is  significant  amount  of 
energy  loss  during  the  aforementioned  equalization  processes.  To 


complete  the  same  trip,  the  average  SOC  drops  by  33.74%  in  Case  1, 
while  in  Case  2,  by  37.16%.  The  corresponding  WOR  is  1.102. 

Simulation  has  also  been  performed  to  the  situation  when 
battery  module  equalization  occurs  not  from  the  start  of  the  trip, 
but  rather  in  the  middle  of  the  trip  as  the  SOC  difference  increases. 
For  this  purpose,  we  consider  the  same  battery  configuration  as  in 
Table  2,  while  all  the  modules  have  equal  starting  SOC,  i.e.  80%.  The 
SOC  deviates  for  each  module  during  the  first  645  s,  and  thenafter 
the  equalization  process  starts.  As  shown  in  Fig.  10,  the  SOC 
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Fig.  14.  Estimation  of  battery  module  internal 
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difference  between  battery  1  and  battery  2  was  2.05%  and  that 
between  battery  2  and  battery  3  was  1.53%,  respectively.  With  the 
equalization  process,  the  modules  are  equalized  to  50.45%  by 
1050  s.  Then  the  whole  battery  pack  completes  the  driving  cycle  at 
SOC  41.95%.  The  total  energy  required  to  travel  the  trip  was  38.05%, 
and  the  WOR  was  1.105  in  this  case. 

The  above  simulation  results  show  that  equalization  can  handle 
the  presence  of  degraded  module  during  vehicle’s  operation 
however  there  is  a  considerable  amount  of  SOC  loss  during  equal¬ 
ization.  Highly  degraded  modules  have  significantly  high  internal 
resistance,  which  would  waste  the  battery  energy  and  in  turn  lead 
to  uneconomical  operation.  If  the  module-specific  internal  resis¬ 
tance  can  be  identified  online,  the  vehicle  owners  can  perform  the 
economical  analysis  based  on  their  preferred  trips,  and  make 
reasonable  decision  on  module  replacement. 

5.3.  Subspace  parameter  estimation 

The  subspace  identification  scheme  described  in  Section  4  is 
applied  to  the  two-module  case  for  estimating  the  battery  internal 
resistances.  As  described  earlier  in  battery  modeling  section,  there 
are  surface  internal  resistance,  load  resistance  and  terminal  resis¬ 
tance  that  appear  in  the  state-space  model  of  battery  equalization. 
This  study  mainly  considers  the  following  parameters:  surface 
capacitance  (Cs),  bulk  capacitance  (Q),  inductors  (Li  and  L2 )  in  the 
equalization  circuit,  capacitors  (Co)  from  the  equalization  circuit, 
surface  internal  resistance  Rs,  and  bulk  internal  resistance  Ri.  The 
nominal  values  for  the  capacitance  and  inductance  parameters  in 
the  two-module  simulation  case  are  set  as:  Q  —  90000  F, 
Cs  =  9000  F,  Li  =  L2  =  13  mH  and  C0  =  12  pF. 

In  order  to  assure  the  unique  solution  or  convergence  of  the 
parameter  estimation,  the  condition  of  persistent  excitation  (PE)  is 
evaluated  for  the  input  signals  [29],  The  PE  order  of  is  n  if 

'  c(0)  c(l)  •••c(n  —  1) I 

c(l)  c(0)  -c(n-2)  (26) 

c(n  -  1)  c(n  -  2)  -c(0) 

is  nonsingular,  where  the  covariance  of  input  is  given  by 

c(k)  =  fim  j  J2  u(i)u(1'  -  k )  (27) 


Determination  of  PE  order  is  based  on  the  evaluation  of  the 
covariance  matrix  and  check  if  it  is  full  rank.  The  number  of 
parameters  to  be  estimated  should  not  higher  than  the  order  of  PE 
that  can  be  achieved  for  the  given  input-output  data.  Using  eigen¬ 
values  is  a  more  reliable  method  to  check  the  matrix  singularity. 
The  condition  number,  i.e.  the  ratio  of  the  largest  eigen-value 
magnitude  to  the  smallest,  is  a  more  trustable  index  to  just  if 
a  matrix  is  full  rank  even  though  all  eigen-values  are  non-zero.  If 
the  condition  number  is  very  large,  the  full  rank  condition  is  not 
practically  valid. 

In  this  study,  the  input  is  the  duty  ratio  of  the  equalization 
circuit.  The  input  data  segment  for  every  100  s  is  used  to  evaluate  if 
the  4th  order  PE  is  valid.  The  condition  numbers  and  the  maximum 
magnitude  eigen-values  shown  in  Fig.  11  validate  the  4th  order  PE 
condition. 

The  internal  resistances  for  Battery  B-l  (Rest,  Rdi  and  Roi)  and 
Battery  B-2  (Rcs2,  Rci2  and  R02)  are  estimated  using  the  subspace 
identification  method  described  earlier.  The  online  estimation  of 
the  internal  resistances  of  B-l  and  B-2  are  shown  in  Fig.  12.  The 
relative  errors  in  estimating  the  internal  resistances  are  plotted  in 
Fig.  13,  which  are  within  18%.  Such  accuracy  is  generally  acceptable 
for  fault  detection  and  WOR  evaluation  purpose.  Notice  that  the 
estimation  error  for  Ro  is  relatively  smaller,  while  that  for  Rcs  is 
relatively  larger. 

Similar  simulation  study  is  then  performed  for  the  3-module 
equalization  process.  In  this  case,  the  parameters  we  identify  are 
Rest.  Rdi  and  R0i  for  Battery  B-l,  Rcs2,  Rci2  and  R02  for  Battery  B-2, 
and  RCS3,  Rd3  and  R03  for  Battery  B-3,  respectively.  The  parameters 
Rcs  and  Rci  are  0.001  Q  and  0.01  Q  whereas  Roi  —  0.054  Q, 
R02  —  0.081  Q  and  R03  —  0.108  fi.  The  three  plots  in  Fig.  14  show  the 
estimation  results  for  the  internal  resistances  for  each  of  the  three 
battery  modules,  respectively.  Fig.  15  shows  the  relative  errors  for 
these  estimations.  For  Battery  B-l,  Rcsi  has  mean  deviation  of  0.8% 
with  variation  from  -7.8—9.25%,  Rcu  has  a  mean  deviation  of -2.5% 
with  a  variation  from  -7.8  to  1%  and  Roi  has  mean  deviation  of  4.5% 
with  variation  from  -2.75—11.75%.  For  Battery  B-2,  Rcs 2  has  mean 
deviation  of  1%  with  variation  from  -4.6-9.4%,  Rci2  has  a  mean 
deviation  of  -2.2%  with  a  variation  from  -9.6  to  2.66  %,  and  R02 
has  mean  deviation  3.9%  with  variation  from  -2.4—11.15%.  For 
Battery  B-3,  Rcs 3  has  mean  deviation  of  1.7%  with  variation  from 
-8.9-9.32%,  Rci3  has  a  mean  deviation  of  -1.15%  with  a  variation 
from  -8.6-8.8%  and  R03  has  mean  deviation  of  1.6%  with 
variation  from  -3.0-10.45%.  For  all  cases,  the  relative  errors  of 
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Fig.  15.  Relative  error  in  estimating  battery  module  internal  resistance  for  the  3-module  case. 
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estimation  are  bounded  by  11.75%,  which  is  again  acceptable  for 
degradation  evaluation  purpose. 

6.  Conclusion 


:s  both  the  ei 
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This  paper  presents  two  ideas  around  I 
dation  of  battery  pack  for  electric  vehii 
resistance  of  each  individual  module  is  e 
vehicle  operation  using  the  battery  eqi 
identification  method  is  adopted  for  sc 
the  module-wise  internal  resistance  r 
specific  WOR  can  be  evaluated,  which  inc 
due  to  high  internal  resistance  of  certair 
performance  recovery  if  a  healthy  module  can  be  ir 
lation  study  has  been  performed  using  an  actual  driving  cycle  for 
a  commuting  trip.  The  identification  of  the  internal  resistances 
demonstrated  less  than  18%  relative  error,  which  is  acceptable 
performance.  As  battery  module  replacement  is  a  big  issue  for  EV 
maintenance,  the  proposed  ideas  provide  a  convenient  way  to 
onboard  health  map  identification,  along  with  a  quantitative 
measure  for  evaluating  the  necessity  of  replacement  for  certain 
module. 


Appendix.  Battery  equalization  model  derivation 

The  schematic  diagram  for  the  battery  equalization  circuit  is 
shown  in  Fig.  2.  Fig.  A.l  shows  the  circuit  diagram  for  battery  1 
when  switch  Qj  is  in  on  state. 

Ro,  R“  R|-1  l' 

r^yn— n 

Ei  CS,  iL,  Cli 

Qi 


Fig.  A.1.  Equivalent  circuit  diagram  for  battery  equalization  when  Q,  is  on. 
The  relevant  dynamics  can  be  modeled  as: 


Vest  —  r  !ti  —  d  r  ^csl 
*-csl  Kcs1  t-csl 

(A.la) 

1  •  1  1/ 
vcLl  =  r — H.1  ~  p  r — *csl 

f-cLl  kcL1<-cL1 

(A.lb) 

in  —  fin  -  Vest  -  vcu 

(A.lc) 

V0  =  ^in 

C-0 

(A.ld) 

The  circuit  diagram  for  battery  1  when  switch  Qi  is  in  off  state  is 
shown  in  Fig.  A.2,  for  which  the  associated  dynamics  can  be 
modeled  as 

_ T<WV-|  |-VjvV-| rH 

Cs,  cti 


Fig.  A.2.  Circuit  diagram  for  equalization  process  when  switch  Qi  is  off. 


^S1  -c^iu~R^Vc* 

(A.2a) 

VcU  =  Ql,  ''L1  ^  Rcl,  Ccl,  Vcs1 

(A.2b) 

in  =  ‘n  -  L]  ucsi  -  L]  vcu 

(A.2c) 

V0  = 

(A.2d) 

If  the  duty  ratio  of  the  circuit  switch  is  denoted  by  d,  an  averaged 
model  can  be  obtained  by  combining  Eqs.  (A.l)  and  (A.2),  i.e. 
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